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Abstract 

We find an effective Hamiltonian describing the process of second-harmonic gener- 
ation in the far-off resonant limit. We show that the dynamics of the fundamental 
mode is governed by a Kerrlike Hamiltonian. Some dynamical consequences are 
examined. 
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1 Introduction 



The superposition principle plays an essential role in understanding the con- 
ceptual foundations of quantum mechanics. Stated in physical terms, it assures 
that we cannot speak of an objective state of a system prior to a measurement. 
In fact, one should admit that the system may be described as a superposi- 
tion of macroscopically (or mesoscopically) distinct states, often known as 
Schrodinger cat states [1]. 

Therefore, the methods for the generation of such superposition states are of 
fundamental interest. Among other systems, vibrating molecules [2] or crys- 
tals [3], trapped ions [4], and Bose condensates [5] have been proposed as 
candidates for doing the task. In fact, mesoscopic Schrodinger cat states for 
the vibrational motion of ion traps have been experimental realized [6]. 
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The issue of generating optical Schrodinger cats (namely, superpositions of 
distinguishable states of the electromagnetic field) has attracted as well a 
great deal of attention [7]. For example, Brune et al [8] have shown that 
the conditional measurement on atoms exiting a high-Q cavity may force the 
radiation inside the cavity to be in a catlike state. The same basic concept 
of conditional measurement has also been suggested for entangled states [9], 
nonlinear birefringence [10], or self-Kerr phase modulation [11]. The state 
collapse due to photodetection has also been suggested as a source of cat 
states [12]. 

The Kerr effect provides a nonlinearity of particular interest for generating 
field cat states [13]. Apart from their intrinsic simplicity, Kerr-based schemes 
have the specific advantage of not relying on conditional measurements. How- 
ever, the realistic values of Kerr coefficients are quite small, thus requiring a 
large interaction length. Then, losses become significant and may destroy the 
very delicate quantum superpositions. In short, although very appealing from 
the physical viewpoint, Kerr schemes are not generally considered to be realis- 
tic. In spite of this belief, there has been a very recent proposal of generating 
a Kerr kitten that serves as a quantum seed for suitable amplification [14], 
that could become fully feasible. 

On the other hand, it should be noted that perhaps the simplest nonlinear 
optical process is second-harmonic generation, which in addition exhibits a 
rich spectrum of nonclassical features. Considerable attention has been paid 
to issues such as photon antibunching, squeezing, or collapses and revivals [15]. 
However, the generation of Schrodinger cat states has gone almost unnoticed. 
Perhaps, a relevant exception is the work of Nikitin and Masalov [16], showing 
that, at resonance, the quantum state of the fundamental mode evolves into 
a superposition of two macroscopically distinguishable states. This point may 
be considered as well rather unrealistic, because the stringent experimental 
constraints of having perfect phase matching and low decoherence could be 
rather difficult to attain. 

In this paper we pursue the dispersive limit of second-harmonic generation, 
which seems to be almost ignored in the literature. At this respect, our main re- 
sult is to show that, for detunings large enough, the dynamics can be described 
in terms of an effective Hamiltonian essentially identical to that governing the 
propagation of the fundamental mode in a Kerr medium. 

Apart from its intrinsic interest, such an identification could be more than 
an academic curiosity: second-harmonic generation is, for a variety of rea- 
sons [17], more robust than the Kerr effect as for noise-limiting factors. In 
consequence, this scheme could be an experimentally feasible proposal to gen- 
erate the yearned optical cat states. 
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2 Dynamics of second-harmonic generation in the dispersive limit 



Second-harmonic generation is usually described (without assuming perfect 
resonance) by the following model Hamiltonian (in units h — 1) 



where a (at) and b (tf) are the annihilation (creation) operators of the fun- 
damental mode of frequency uo a and of the second-harmonic mode of fre- 
quency u b , respectively. When perfect matching conditions occurs, they satisfy 
uj b = 2u a . The constant g, describing phenomenologically the coupling between 
both modes, can be always chosen as real. 

Classically the above problem admits an exact solution, leading to the pos- 
sibility of the energy being transferred completely into the second-harmonic 
mode [18]. Unfortunately, the corresponding quantum dynamics is a touchy 
business, though different algebraic, semiclassical, and numerical approaches 
have been developed (for pertinent literature on these approximations see 



H — if o + H i nt = u) a a 



] a + u b b ] b + g(a 2 b ] + a^b), 



(1) 



Ref. [19]). 



By noting that the Hamiltonian (1) admits the constant of motion 



iV = a) a + 26 f 6, 



(2) 



it can be recast in the following simple form 



a 



N, 



o 



3 



(3) 



int 




where 



A = co b - 2u a 



(4) 



is the detuning. 



It is straightforward to check that 



[Hq, Hi n t] = 0, 



(5) 
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so both Ho and ifi nt are constants of motion. Since Hq determines the total 
energy stored in both modes, which is conserved, we can factor out exp(—iH t) 
from the evolution operator and drop it altogether. 

In order to examine the quantum evolution we shall use a standard numerical 
approach developed by Tanas and coworkers [20] starting from the original one 
of Walls and Barakat [21]. First, we consider the case of perfect resonance and 
we take the typical experimental conditions; namely, the fundamental mode 
in a coherent state with an average number of photons n a and the harmonic 
mode in the vacuum. 

For definiteness, we shall analyze the quantum state by computing the Q func- 
tion. In figure 1 we have plotted a gray-scale contour map of this distribution 
Q(a) for the fundamental mode with n a = 10, and for different values of 
the dimensionless time r = gt^/2n a . Results are in complete agreement with 
the previous work by Nikitin and Masalov [16]: at the beginning of the inter- 
action (r < 2) the uncertainty contour turns to an ellipse, whose long axis 
is perpendicular to the field vector (so the state is squeezed in amplitude). 
Further interaction (r ~ 4) results in a bifurcation of the Q distribution and 
the formation of a superposition of two macroscopically distinguishable states. 
Subsequent development (r > 5) destroys this state as well. 

Next, we consider the case when some mismatch is present. The previous 
numerical procedure is still valid, but one must take care of the diagonal 
extra terms introduced by the detuning. To check how the detuning affects 
the cats, in figure 2 we have plotted the contours of the Q function at r = 4, 
when the cat should appear, for different values of the dimensionless detuning 
A/g. One easily notices that the cat becomes more and more blurred as the 
detuning increases, until A/g is so large and the process is so inefficient that 
the fundamental mode does not feel any appreciable interaction and the cat 
is lost. Moreover, we have found no numerical evidence of any further cat for 
other times around r = 4. 

This is a reasonable result from the physical viewpoint, since one could expect 
some kind of continuity. However, it differs from the numerical computations 
of Nikitin and Masalov [16], who find that the cat disappears as soon as some 
detuning is present, irrespectively the smaller this detuning is. The reason is 
the way they deal with the off-resonant situation: following a classical argu- 
ment, they assume that the presence of detuning translates into the fact that 
the coupling constant has the time dependence g = \g\ exp[(n a — nb)ut], where 
n a and rib are the refractive indices of the medium at the frequencies oo a and 
ujb- In such a way, the interaction Hamiltonian H[ nt becomes time dependent 
and the numerical method becomes rather complicated, at difference of what 
happens with our very simple model. 
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From now on, we shall be interested in the dispersive limit of the model gov- 
erned by the Hamiltonian (1). In other words, we shall consider the situation 
when 

|A| >#(n a + l){n h + 1), (6) 



where n a and n& denote the average photon numbers in the modes a and 6, 
respectively. It is worth emphasizing that this is a limit physically realizable 
in practice [22]. Now, the essential point is that g/A can be considered as a 
small parameter and, according to the method developed in Ref. [23], we can 
apply a unitary transformation to the interaction Hamiltonian (3) 



H eS = UH int U\ 



(7) 



where 



U = exp 



9 



(a 2 6t - a t 2 6 ) 



(8) 



To gain a deeper insight into the physical ideas underlying the method we can 
introduce the operators 



X + = fcV, X_ = a f2 6, 



X 3 = -(tfb-a!a), 
that satisfy the commutation relations 



(9) 



[X 3 ,X ± ] = ±X ± , 
[X + ,X_]=P(X 3 ), 



(10) 



where P refers to a quadratic polynomial function of the diagonal operator 
X 3 . Then, while the first relation is the same as the corresponding one for the 
standard su(2) algebra, the second one defines what is called a nonlinear or 
polynomial deformation of su(2) [23]. From this viewpoint, the transforma- 
tion (8) can be recast as 



U = exp 



9 



{X+-X-) 



(11) 



For su(2) operators this would be simply a rotation. Here, it represents a 
small nonlinear rotation with the parameters chosen in such a way so as to 



5 



cancel [up to order (g/A) 2 ] the nondiagonal terms appearing in the trans- 
formed Hamiltonian. Moreover, its unitary character ensures that H e g has the 
same spectrum as H mt (so they are physically equivalent). In consequence, by 
expanding Eq. (7) in a power series and keeping terms up to the order (g/A) 2 , 
we finally get 



H cS = -(6+6 - a+o) - A [46+60+0 - (a+a) 2 



(12) 



with 



\- 9 
A' 



(13) 



The effective Hamiltonian (12) describes the dispersive evolution of the fields, 
but the essential point is that it is diagonal, which implies that there is no 
population transfer between the modes (as it is expected in the far-off resonant 
limit). The first term in Eq. (12) does not affect the dynamics and just leads 
to a rapid oscillation of the wavefunction. It is worth noting that the small 
rotation generated by the transformation (8) just adds small corrections of 
order g/A to the physical observables and, thus, the state vector need not to 
be transformed. 

Now, it is easy to find evolution of the density matrix of the fundamental 
mode. When the harmonic mode is initially in the vacuum, the two terms 
containing b^b in Eq. (12) do not contribute. In addition, the linear term in 
a+a lead just to a c-number phase shift and can be also omitted. Then, the 
effective Hamiltonian reduces to 

H K = A(a+a) 2 , (14) 



which is nothing but the interaction Hamiltonian that governs the state evo- 
lution of the single-mode field a in a Kerr medium. Then, the evolution is 
factorized all the times (in the frame of this approximation) and, in conse- 
quence, starting from a coherent state in mode a, one will recover M copies 
of this initial state at times At = n/M [13]. In particular, for At = n/2 the 
density matrix of the fundamental mode becomes 

PaW2\) = \ll> c ) aa (ll>c\, ( + 5) 



where \ip c )a is a superposition of two coherent states with opposite phases; 
i.e., 



V2 



3 i7r/4 



a) a + e 



-i7r/4 1 



a) 



(16) 
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To check these predictions we have numerically calculated the Q function from 
the exact Hamiltonian (3) at the time t = 7r/(2A). The initial coherent state 
in mode a has n a = 10 and we have assumed a detuning A/g = 50, which 
is enough to guarantee the validity of our approximation (6). In figure 3 we 
have plotted the results of the computation, showing an apparent excellent 
agreement with the expected cat. However, the Q function by itself is not 
very suitable to distinguish between coherent and incoherent superpositions. 
To further confirm the coherent behavior we have considered the fidelity 

Ht) = a(4>MtM c ) a (17) 



between the ideal two-component cat density matrix and the state coming 
from realistic evolution. In figure 4 we have plotted the resulting fidelity as 
a function of the rescaled time gt for A/g = 50. We can see the presence of 
small peaks corresponding to the presence of M— component catlike states, as 
well as the strong peak at the time At = %/2 predicted by our theory when 
M = 2. 

When the initial state of the harmonic mode is in a superposition of num- 
ber states instead of the vacuum, the dynamics of the fundamental mode is 
drastically affected. For simplicity, we shall take an initial coherent state in 
both modes, namely, \a) a <8> \P)b- Then, we obtain for the density matrix of 
the fundamental mode the following expression 



oo 

p a (t) = e-H 2 ST e ^(n 2 -m' 2 ) exp [|0|2 ^ _ e -4iXt(n-m) 



a n a* m 



x r -— \n) aa (m\. (18) 
V n\m\ 

The entanglement between modes a and b is evident, as well as the fact that 
only when At = n/2 they become disentangled and two copies of the initial 
coherent state in the fundamental mode will be created. Indeed, the density 
matrix of the total system at such time is 

p(Tr/2\) = \(3) bb ((3\®\^ c ) aa (^ c \, (19) 



and only a two-component cat can be produced in this case, although with 
lesser fidelity than for the vacuum. On the contrary, when At = 7r, the initial 
coherent state is reconstructed. 

Finally, some comments concerning the physical feasibility of the proposed 
scheme seem in order. Of course, a simple comparison of the parameters in 
figures 1 and 3 clearly shows that the times required for the appearance of 
the cat in the dispersive limit are an order of magnitude larger than in the 
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usual case of no detuning. The advantage of this dispersive limit is twofold: 
one does not need to take care about phase-matching considerations and, more 
important, if the process takes place in a cavity, one could expect from previous 
results, e.g. in cavity quantum electrodynamics [22], that the quantum features 
of the state may be undestroyed for large interaction times. In addition, the 
width of the peak appearing in the fidelity is of the order of 5 in units of gt, 
while the cat appearing in figure 1 disappears in much shorter times. Thus, 
the required time resolution in this dispersive limit is lesser than the one 
needed for the perfect frequency-matching case. In spite of these advantages, 
the required times are still so long that other decoherence factors could affect 
the final results. 

The effective Hamiltonian (12) has the virtue of simplifying calculations that 
otherwise would be rather involved. To cite only an example, the squeezing 
properties of the fundamental mode are also affected by the initial state of 
the harmonic mode. For initial coherent states in both modes \a) a <E> \/3)b (for 
simplicity, we take a and (3 as real numbers) we easily get for the variance 
a x = ( x2 ) ~ ( x ) 2 °f ^e quadrature x = (a + a})j\[2 



al(t) = l - + a 2 (e- 8z2T2 - e~^ T ") cosfiT 

-2a 2 T 2 (4e- 8z2T2 - e" 4 ^) cos fiT 
-2a 2 T(2e~ 8z2T2 -e- 4z2T2 )sinnT 

+ a 2 (l-e' 4z2T2 )cosQT, (20) 

where T = At < 1, z 2 = a 2 + 4(3 2 , and Q = 4a 2 - 8(3 2 . Then, for sufficiently 
large values of f3 2 squeezing in the fundamental mode disappears, according 
to previous numerical studies [16]. 

One could be tempted to extend the method to other nonlinear optical process 
in order to get effective Hamiltonians that are diagonal in the desired basis. 
In fact, the program can be carried out. For example, for third-harmonic 
generation one obtains, after repeating the same steps as before, 



»cff 



A {96+6 [(a f a) 2 + a f a] - (a f a) 3 - 6(a f a) 2 } , (21) 



but now the final result contains a cubic nonlinearity and the problem is not 
facilitated at all. 
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3 Concluding remarks 



What we expect to have accomplished in this paper is to present a simple 
model for the description of the dispersive limit of second-harmonic generation 
in terms of an effective Hamiltonian that is diagonal and contains Kerrlike 
nonlinearities. Of course, in the practice other off-resonant processes could also 
contribute to the third-order nonlinearity. However, our effective-Hamiltonian 
approach can explain (al least qualitatively) some unexpected and fascinating 
behaviors of the so-far unnoticed dispersive limit. 

Of course, it is well known that when some interaction energy is off-resonant 
then, in a higher order of perturbation theory an effective interaction energy 
can be constructed that gives rise to resonant processes. At this respect, in 
Ref. [23] it has been shown how the method of small rotations can be used to 
tailor resonant interactions from nonresonant processes, just as demonstrated 
in this paper for the relevant case of the second-harmonic generation. 

Finally, it is worth noting as well that the Hamiltonian (12) is similar to that 
obtained in Ref. [24] for the Dicke model in the large detuning limit (when 
the effective Hamiltonian was proportional to S%). 
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Fig. 1. Gray- level contour plot of the quasiprobability distribution Q(a) of the 
fundamental mode (with perfect frequency matching) at different times r = gt\J2n a . 
Initially, the fundamental mode is in a coherent state with n a = 10 and the harmonic 
mode is in the vacuum. In all the cases in the horizontal axis we plot Re a and in 
the vertical we plot Im a. 



Fig. 2. Gray-level contour plot of the quasiprobability distribution Q(a) of the 
fundamental mode for the same initial state as in figure 1 calculated for r = 4 at 
different values of A/g. 



Fig. 3. Q function and the corresponding contour plot of the fundamental mode 
(initially, the fundamental mode is in a coherent state with n a = 10 and the har- 
monic mode is in the vacuum) with a detuning of A/g = 50. The function has been 
calculated calculated from the exact interaction Hamiltonian (3) at Xt = tt/2. 



Fig. 4. Fidelity of the fundamental mode respect to the ideal cat state as a function 
of the rescaled time gt. The state evolution has been calculated from the exact 
interaction Hamiltonian (3). The initial state is as in figure 3 with the same detuning 
A/g = 50. The main peak appears at the time Xt = n/2 predicted by the theory. 
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